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I ' Abstract 

gj ■ Recently we have discussed the generalized parametrized Klein-Gordon 

(3JT)' equation for curved spacetime. We have also discussed its derivation from 

jyj . several approaches, the direct Feynman parametrization, the state func- 

^ ' tion entropy or equivalently the information theory approach, and the 

' rj!, t stochastic differential equation approach. We have even suggested a gen- 

^- . eralization of the statistics of the entropy to the generalized entropies and 

(**] ' derived the particular nonextensive statistics parametrized Klein-Gordon 

O I equation, and discussed its nonlinear FPE replacement of the complicated 

Gibbs-Boltzmann statistics entropy derived analog with complicated non- 
linear potential or drift and diffusion coefficients. In this article we apply 
these previously derived results to the quantum transport in abruptly 
^vj ■ coupled curved space-time heterostructures, applied here specifically to 

Black-Hole event horizon coupling to normal curved space-time. We de- 
rive the coupling self energy, and the Garcia-Molliner surface Green's 
functions from which we can calculate the surface area and entropy. We 
then derive the nonequilibrium transport equations for the radiation from 
p~Z • the Black-Hole. We discuss the theory application to Worm Holes and 

(— "•; ' quantum analogues. 
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1 Introduction 

X 

H ! 2 derivation 

C3 



The parametrized Klein-Gordon equation can be obtained from the max- 
imum entropy equation as the maximization of the entropy < S >= 
— J PlnP with the moments < x 2 >= J x 2 P(x,X)dx multiplied by the 
Lagrange multipliers as 
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S[<S(\)>]+S[J2fr& 2 )]=0 (i) 



The maximization allows us to derive the least biased probability distri- 
bution function 

-X>(*,) 2 

P(x , A) = JV(A)e '=1 (2) 

The Lagrange multipliers are evaluated by the identity —-^§§- —< (Xi) > = 
2Di\ . The partition function is the inverse of the normalization Z — i. 
The PDE of the evolution of the least biased distribution Eq.@ is for the 
parametrization — /? < A < /? 

gfjf dp 0i P ,, 

d\ 2^ dx t 2 ' [ ' 

The Lagrange multipliers which are the inverses of the diffusion coefficients 
are when set to the coordinate Di=i,2,3 = 2, D4 = ■% and upon taking the 
Laplace transform (or Fourier transform with A— > — i\ with appropriate 
analytic continuation) P(x, A) = J P(x, s)e~ Xs ds becomes with s = — ' £ 
the massive Klein-Gordon equation 

dP_ dP_ dP_ 1 dP _ m 2 c A 

dx^ + d^ + d^ + ^ d{-itf ~ ~^~ () 

The generalization to the curved space time is done by the inclusion of dif- 
fusion coefficients from the metric tensor. The diffusion coefficients are the 
inverses of the diagonal metric tensor. Alternatively the curved spacetime 
potential is utilized, the Shroedinger-like PDE is obtained, a transforma- 
tion to the Fokker-Planck PDE is made and therefore drift and diffusion 
coefficients are obtained. Utilizing the Girsanov theorem transforms that 
in turn to a drift free form which has the coordinates transformed as the 
drift dependent coordinates Alternatively SDE transformations can be 
made. 

The diffusion coefficients we will be working with are the diffusion 
coefficients from the inverses of metric tensor of the Black-Hole curved 
space-time. 

The Schwarzchild metric in Cartesian coordinates is 

dT " = - n+ffj^ 2 + C 1 + ^7) 2 ( d ^ + d V 2 + dz ") (5) 

The metric space- like differential equation is , with t* — —it 

dr = A x (x)dx 1 + A y (x)dy j + A z (x)dz k + A t (x)dt* i (6) 

with the stochastic differential equations 

dr = dW{\) (7) 
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This equation when the diffusion coefficients are multiplied to the left and 
the metric product taken becomes the space-like metric, and the delta 
correlated noise terms (dW(\)) 2 — d\ of the space-like metric vector dr. 
The delta correlated noise terms (dW t {\)) 2 + (dW x (\)) 2 + {dW y {\)) 2 + 
(dW z (X)) 2 = {dW(\)) 2 = d\. That is the sum of delta correlated Wiener 
processes is a Wiener process, and the square of the delta correlated pro- 
cess is replaced by its average d\ [2]. 

The SDEs are utilized to obtain the Fokker-Planck equation [2] for the 
process, which is 

9P&V _ £ A _9% MG a D t d 2 (1+^f) 2 

0A 2 ^ fla* 3 U c 2 r ; jjr+ 2 a(-Ji) 2l (l-||2)2^ 

(8) 
Therefore the diffusion coefficients, with the D — 2, D t — -% factored 
out from the noise terms and written explicitly. The constant diffusion 
coefficients are now generalized to the inverses of the Schwarzchild metric. 
The Fokker Planck equation is solved by several methods, and we write 
the short time transition probability solution [3] 

The solution of the time here evolution parameter A dependent Fokker- 
Planck equation is obtained as 

-(x-x')[rT^ 1 (x' ,t')](x-x') 

P(x, t; x', t') = (9) 

The FPE equation is also to be transformed by SDE transformations 
and compared to the solution Eq.((9|. We write the SDE for example the 
X coordinate as 



We set the diffusion coefficient of this transformed SDE to a constant and 
obtain after a change of variables x' — r'cos(0')sin(</>') the differential 



X = D " I \/(l-S?) 2 ^'= (ID 
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The SDE can be written as 

dj y * I ^7> {1 _^y2 dx + ^ m jJ dw ^ X) (12) 

Z ox c A r 

and the transformation of coordinates for dy,dz,d,T is made similarly to 
obtain the drift coefficients 

° (1 _M& ) (l i MG\2 
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This stochastic process has the Fokker-Planck evolution, and we redefine 
£* = — it => — zc£ and x— > x, y, z, t* 

T->4t£ 

The two point transitional probability for short times solves the Fokker- 
Planck equation with the transition time as the difference A — A' = AA 

(i-£'-l(i , l t")A)| i 

P(x,t*;x,t*') = - //"?" —• ( 15 ) 

3 Derivation from potential 

the equations for the Schwarzchild metric and their solution have been 
discussed in terms of the derived drift and diffusion coefficients. An al- 
ternative method is to derive these equations from the state function en- 
tropic measure, with the gravitational potential due to the Black-Hole a 
constraint . We derive the least biased probability distribution for the 
potential < V >. 

The potential is obtained from a Lagrangian L — T — V , with the 
variation obtaining the Euler-Lagrange equation -§^t^ — §b- = 0. The 

dr 

differential for example %£. = ^%r = -. — m£, 2 W? an d the differential 
The potential then for the example x coordinate 



(1 i wiy ^ 

_d_de_ __ { - L+ ^7> d 

df dr (1 M .fy. dt' 

becomes £ro£. 

The kinetic energy has a mass term, and as we parametrized the mas- 
sive Klein-gordon equation, we will need to rewrite the mass term as the 
parametrization A. ...and note that other parametrizations could have been 
made which after calculation could be taken as a limit approaching zero, 



and as the massive Klein-Gordon equation utilized in the photon propa- 
gator was parametrized and after calculation the mass term was taken as 
a limit approaching zero. 

The kinetic energy term T = ^(^) 2 has a mass term, which for exam- 
ple can be written as T = t"\// -§rze~ sX d\. The Euler-Lagrange equa- 
tion has the force relationship m4-(^r) = 4p- which is obtained from the 
discussed metric and from which we derive the relations Jj [- — 4&T2 (at* 7 — Mfl <i ^+ 



MO,2 V x + n- MO >2 a 'A — F* 



L(i_ jw<i)2 \at* (1 a )" 

(i+«£) 2 

n — A f , u 1 2 " x ' f 1 ^ c t T ^ 

The very complex force components are simpler in the spherical polar 
coordinates due to the radial dependence of the potential. We write the 
potential as the radial V = G , M . We can rewrite this in Cartesian 

coordinates. The velocity squared is v 2 = v x 2 + v y 2 + v- 2 + v t * 2 , and 
the radial coordinate is the modulus of the Cartesian two body distance 
r = y/(x - x o y + (y- y ) 2 + (z - z ) 2 + (t* - i* D ) 2 . 

The entropy is the, not to be confused with the spacelike metric ds 2 , 
< S(\) >= - / PlnP, the potential < V(x, A) >= / V(x)P(x, \)dx. The 
maximization yields 

- J B[(af_af)) 2 + ot y(a',A')] 

P(?,X;xX>) = _ (16) 

This can be rewritten by uncompleting the square of (x— x — a(x , A', \)a(x , A')), 

— fl[(x — x —cn(x ,X ,X)a(x ,A ))] 

p(*, a; *V) = ^ — (it) 

We compare this to the short time transition probability solution Eq. (|15p 
and note that the Lagrange multipliers are P = 2 da\ , the partition 
function Z = \ALrDAA and the Lagrange multiplier we factored out, 
&{x' , A', A) a complicated function we can obtain from the Fokker-Planck 
equation or by direct integration, and by the identity — i§^- =< V >. 
Also we note that the flat spacetime operators we utilized can be modified 
to include curvature by the similar change in coordinates as in our previous 
section's derivation. 

Another approach is to rewrite the moments of the coordinates as the 
Hamiltonian operators p = -J^z + V 2 , and obtain 

P(x,(3) = e-^ + * 2+aV ( 3 '' X ' )] (18) 

The partial differential w.r.t. A obtains for — j3 < A < /3 

i — 1 

The Fourier transform i J"] P(x, iuj n ) e „" and J P(k,icj n ) trans- 


forms the Shroedinger-like PDE when in the two-point form with a delta 



function source to the Matsubara Green's function 

P(k, iw„) = — ^ (20) 

icu„ — e (k) — aS(k, iu>„) 

The inverse Fourier transform, with the potential self energy equal to zero 

(x-x') 2 

for simplicity, is the distribution function P {x, A; a?', A') = e /4 2 "^, A A A ■ 

The potential energy can be included via the transformation of coor- 
dinates K 2 = k 2 + aV(k,icu n — 0) ...the inverse Fourier transform yields 

(X-X 1 ) 2 

a Gaussian P(X, A; A"', A') = e . , and the coordinate transforma- 

V4ttDAA 

tion is then to transform X — > x by direct Fourier transforms, an al- 
ternative to PDF transformations, coordinate transformations, and SDE 
stochastic transformations. 

A point of interesting discussion here is of the discrete Fourier trans- 
form and the Matsubara frequencies. The sum corresponds to even fre- 
quencies or Bosonic Bose-Einstein statistics, and odd frequencies or Fermionic 

Fcrmi-Dirac statistics. In this case the f(k,iuj n — > co+iu) — — - — 4- — ■ — The 

point here is that the momentum vector is a 4Dimensional relativistic vec- 
tor that includes the real time variable. It is derived from a boson equation 
that is parametrized with an evolution parameter A, here related to the 
physical temperature by the Matsubara derivation. The generalized 5Di- 
mensional equation of 4D+lp lparameter yielding a higher dimensional 
Fermi-dirac statistics and higher dimensional Bose-Einstein statistics is a 
novel feature that will be investigated further in future work. 

The derivation of this curvature potential acting upon Bosonic parti- 
cles here photons has been instructive in pointing out the connexions to 
our entropy state function derivation, the stochastic differential deriva- 
tion, and to this 5D Shroedinger-like relativistic quantum mechanics. 

The real curved space-time equation in our derivation would be built 
from the flat spacetime momentum and therefore wave equation deLam- 
bertian, then modified to include the curvature terms as diffusion coef- 
ficients obtained from the curved space-time metric. A change of co- 
ordinates as before transforms Xw=i Bi(x, A)V 2 ; S + B t * (x, A) V 2 t* — > 
Xw = i ^ 2{ - + V 2 t* here we write the nabla as a partial derivative for 
simplicity of notation. The operator form of the momentum is utilized 
in our maximum entropy derivation, the Matsubara Fokker-Planck or 
Schroedinger-like 5D or 4D+lp equation obtained, and Fourier trans- 
formed to the momentum-frequency space. The subsequent inverse Fourier 
transformation obtains statistics of Fermions and Bosons and if need be 
mixed and fractional statistics. We then have side-stepped the Dirac prob- 
lem of quantizing curved space-time operators. Commutation relations 
are not considered here as problematic as the transformations of curva- 
ture terms which modify flat space-time operators transforms the problem 
entirely to an effective linear flat space-time equations with nonlinear po- 
tential and diffusion transformed drift potential terms. A new flat-like 
space-time which is effectively quantizable, effectively solvable in analytic 
form. Once the quantum mechanics albeit in 5D or 4D+lp spacetime is 
obtained from the flat spacetime operator representation, the curvature is 



restored by the analytic forms of the coordinate transformations which it 
has been discussed in detail can be obtained from the state function, the 
macroscopic PDF or the microscopic SDE levels. The remaining problem 
of interpretation of the parametrization, the evolution parameter lambda, 
and the Lagrange multiplier /3 is the decades old noted correspondence be- 
tween variance, temperature, diffusion and with properly applied analytic 
continuation, quantum mechanics. The theory then becomes an entirely 
well founded physical approach to the quantum statistical mechanics of 
curved space-time in terms of relativistic thermodynamics. We need not 
discuss in detail the correspondence between the maximum entropy prin- 
ciple, information theory and thermodynamical variational theory. 

4 Exact coupling of Normal Curved Space- 
Time and Black-Hole Interior. 

The Hawking problem in obtaining albeit successfully the radiation evap- 
oration dynamics of Black-Holes was to match boundary conditions be- 
tween the normal curved space-time and the Black-Hole interior. The 
problem was complicated by a noted boundary conditions matching pathol- 
ogy in elliptical PDEs. A problem or pathology that does not exist in our 
5D or 4D+lp Schroedinger-like relativistic mechanics of the parabolic 
PDE form. The theory of boundary conditions matching has been dis- 
cussed in detail by T.E. Feuchtwang for the 3D Cartesian and spherical 
polar and generalized coordinate systems. We adhere to this approach in 
our +1 dimensional theory, and note that the (curvedSpacetime)-(Black- 
holclnterior) interface heterostructure is a simple two component device in 
this theory. More complicated strucures would be (infiniteLead)-(Black- 
Holelnterior)-(infiniteLead') and here 'lead' is spacetime, these more com- 
plicated heterostructures of time bias or space bias driving transport of 
translation in time-space.... note that this more complicated structure is 
called a worm hole and the time and/or space bias driving transport of 
translation in time-space is of the nonequilibrium transport for large scale 
phenomenon and quantum mechanical nonequilibrium quantum transport 
for mesoscopic and nanometer and sub nanometer scales.. At quantum 
scales these structures are quantum black holes and quantum worm holes, 
and at Planck scales this more complicated heterostructure is called a 
GEON after the work of John Archibald Wheeler, for which if a metastable 
heterostructure of statistical Boson or Fermion or mixed statistics can be 
analytically derived a theory of fundamental quantum particles can be 
rigorously obtained. 

The simple heterostructure of normal curved spacetime in abrupt cou- 
pling to a Black-Hole is obtained by the delta function coupling method. 
This is obtained from a Green's theorem derivation as 

inf „ 

y / g(x; x"; iuj„)[ihuj n + Vj» + H'(x")]G(x"; x ; iui n )dx" (21) 



inf 

■E 



' G(x; x"; iu n )[ihuj n + V%„ + H'(x")]g(x"; x; iu}„)di 



These equations can be integrated by parts to obtain the Dirichlet 
g(x)\g = toundary — and the Von Neumann boundary condition Green's 
functions y_f' \s-boundary — 0. Other mixed boundary conditions are 
also possible. The derivation is explicitly discussed in [3]. The coupling 
due to matching boundary conditions can be generalized to delta function 
coupling Hamiltonian as h(x) = [-§=,5{x — B — 77)]+ — [-^,8(x + B — r/)] + . 
The curved spherical Black-Hole encourages us to work in spherical 
polar coordinates. We will do that in subsequent derivations, we wish now 
to point out that the surface Green's functions are immediately obtained, 
they are for Dirichlet boundary conditions as an example [JJ 

r , , M _ ( g st (r,r',iuj n )g bh (r,r',iuj n ) 

[ ' >"*>"*=*=** - \[g at {r,r',iw n )+g bh {r,r',iu n )] 

(22) 
and a similar expression in terms of the 1st order partial derivatives 
is obtained for the Von Neumann boundary conditions... these surface 
Green's functions are called the Garcia-Molliner Green's functions. The 
overall Green's function is written in terms of the normal curved spacetime 
Green's function g a t and the Black-Hole Green's function g b h which both 
Green's functions exist throughout space-time yet are restricted by the 
boundary matching we have discussed to the regions of the r 8t > R b h > 

Tbh > + 77 

The Green's function is the probability distribution function for ana- 
lytically continued quantum mechanical functions, and is the probability 
function for hermitian conjugated functions when complex valued. The 
Matsubara Green's function is a real valued function as can be noticed 
from the 5D Schroedinger-like PDE, the solution of which is the Gaussian 
we have discussed before, the entropy for the Black-Hole is calculated 
as < S(/3) >= J G(r,r')lnG(r,r')drdr' . The temperature is set after 
the inverse Fourier transform as discussed in the previous section, and 
with choice of Bose-Einstein , Fermi-Dirac or fractional Haldane statistics 
[6]. The surface entropy is obtained from the surface Green's function as 
S(0)^GlnG\ r=r/=Rbh . 

The overall Green's function can be written as 

G(r,r\ic,,0=5^r\^„)-^-[ gst (r,r'')^ 7 G(r'',r')|^:^^ (23) 

The Green's function can be written as a function of the coupling 
self energy. The Green's function with the transpose equation can be 

evaluated as dl,,,,G\ — —i—-, — ,„ „, } ,„ „., . The Green's function 

can also be written as a function of the coupling self energy, where the self 

energy is E r = — ^ (70a — "fa) r ,g B) ■ The self energy is dependent on 

the method of partitioning the Hamiltonian, and the Dirichlet boundary 
conditions reduce this to [T] YZ. = -^- = 4r^, and where a. = st,bh 

and k a — yhijj — Ti{k a ,iu) — E^ — it], and here we have written it in the 
analytically continued form ioj n — > u + irj. 



A possible simplification can be made that focuses on a region of space- 
time, say the curved normal spacetime, perturbed by coupling to the black 
hole. This is a transformation that can be derived from the atomistic 
discrete matrix method, and as we have shown is also a transformation 
that can be derived from the continuous method. The Green's functions 
are then the Green's functions derived perturbed by the coupling self- 
energy such that Gst = — jn — ^rr ; — ^ ^ — and similarly for the 

Green's function of the black hole. 



5 Nonequilibrium transport 

The equations for the nonequilibrium Green's functions are 

G + = (l + EG r )< ? + (l-E a G a ) + G r E+G a (24) 

The current is a characteristic of the nonequilibrium transport. In 
quantum mesoscopic, nanometer scale interface devices, the device is char- 
acterized by the I-V current-voltage. The device current flow is driven 
by bias, by current, voltage bias, material chemical potential differences, 
here this is a space-time derivation. What would the driving correspond 
to? The driving is a mismatch in potential energy between the curved 
space-time and the black-hole. The analogy to current driving here is 
a current of quantum particles incident on the blackhole...this can be a 
current of virtual particles. The analogy to chemical potential is also 
available for us to utilize in this derivation of the quantum current. This 
can be seen as the simple inclusion of the number of particles as a non- 
constant and controlled by the chemical potential as a constant as in 
< n > a = hik,,-.. — rr- and as derived from Grand canonical ensemble 
considerations. Another intriguing analogy is that here these Fermi-Dirac 
and Bose Einstein statistics are 4D+lp and obtained from the analyt- 
ical continuation of the Matsubara frequencies summation to the inte- 
rn/ 
gration of frequency i ^ g(iu) n )— > J f(uj)g(uj)^ and f(oj) is Fermi- 

71 = 

Dirac or Bose-Einstein as the even or odd frequency summation [7|. This 
can be generalized as by Datta and Maclleland \E\ , as the nonequilib- 
rium Green's function can be written as g + {k, to) = —2if(hio — fi)a(k,uj) 
for the unperturbed Green's function, and as G + (x,x',uj) = —2iF(huj — 
fl(x, x 1 , Ttui) )A(x, x 1 , Tuo). Here the generalized chemical potential is 4vec- 
tor dependent, and the driving bias due to the time and space curvature 
can be seen clearly. This is perhaps a more natural representation, as for 
more complicated structures, say worm holes, the space can be the same 
(?) however the time can be different, with the bias or chemical potential 
being different due to the time difference between the two sides of the 
worm hole. 

The current is obtained from the nonequilibrium Green's function as 

< J(r) > \ r=B = J (I - | 7 )G+(r,r',u,)|w->^ (25) 
= f d r G r (r, S)| r=s E + G°(S, B) - G r Z + d r ,G a (B, r')|r'= S fj 



and here we set the the radial coordinate to the event horizon black hole 
radius for definiteness. 

We see then that the current is nonzero from the following discus- 
sion. The NEGF in the approximation G + — G r S + G a when differenti- 
ated as in the representation Eq. (|23p and its transpose derives the sec- 
ond order partial mixed differential d^ r , =B G r ' a ' + '~ — const / [g st (B, B) + 
g bh (B, B)] r ' a '+>-, and d r , =B G{B, r') = g(B, B)d 2 r „ y G(r", r')\ r „ <r , =B and 
as the Dirichlet boundary conditions rid the uncoupled Green's functions 
d r / r i = B ga(r, r') = 0, the nonequilibrium green's function becomes a sum 
of uncoupled Green's functions only, evaluated at the event horizon radius 
of the black hole. We make the correspondence that this nonvanishing cur- 
rent evaluated at the radius of the black hole corresponds to the Hawking 
radiation. 

6 conclusion 

We have derived an exact transport theory for curved spacetime and 
Black hole heterostructures. This theory is applicable to more com- 
plex structures such as worm holes. We have derived a quantum ther- 
mal Green's functions representation of curved spacetime and utilized the 
4D+lp Schroedinger-like PDEs as the equations of evolution for the rigor- 
ous Green's theorem coupling of normal curved spacetime and the Black 
hole structure interface. We have applied the continuous theory of self 
energy of coupling to the interface coupling of the spacetime and black 
hole horizon and obtained the current evaluated at the radius of the black 
hole. We have shown that even in the approximation of G + = G" E + G a 
corresponding to no large nonequilibrium term contributions the current 
at the black hole horizon is nonzero. We make the correspondence that 
this nonvanishing current corresponds to the Hawking radiation. In the 
future we will derive the full nonequilibrium contribution, thereby the 
change in rate of evaporation of the black hole. Also, numerical simula- 
tions will help visualize the dynamics of this theory, also of interest and 
for future work is the application of this theory to the worm hole physics, 
and the exploration of the limitations if any of the quantum length scales 
of application of this approach to highly curved spacetime. 
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